Voyage to Alpha Centauri: Entanglement degradation of cavity modes due to motion 
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We propose a scheme to investigate whether non-uniform motion degrades entanglement of a 
relativistic quantum field that is localised both in space and in time. For a Dirichlet scalar field 
in a cavity in Minkowski space, in small but freely-adjustable acceleration of finite but arbitrarily 
long duration, degradation of observable magnitude occurs for massless transverse quanta of op- 
tical wavelength at Earth gravity acceleration and for kaon mass quanta already at microgravity 
acceleration. We outline a space-based experiment for observing the effect and its gravitational 
counterpart. 

PACS numbers: 03.67.Mn, 03.65.Yz, 04.62.+V 



I. INTRODUCTION 

A common way to implement quantum information 
tasks involves storing information in cavity field modes. 
How the motion of the cavities affects the stored infor- 
mation is a question that could be of practical relevance 
in space-based experiments 

In this paper we analyse the degradation of initially 
maximal entanglement between scalar field modes in two 
Dirichlet cavities in Minkowski space, one inertial and 
the other undergoing motion that need not be station- 
ary. Our analysis combines, to our knowledge for the 
first time in relativistic quantum information theory, the 
explicit confinement of a quantum field to a finite size 
cavity and a freely adjustable time-dependence of the 
cavity's acceleration. This allows observers within the 
cavities to implement quantum information protocols in 
a way that is localised both in space and in time 0. 
In particular, our system-environment split is manifestly 
causal and invokes no horizons or other teleological no- 
tions that would assume acceleration to persist into the 
asymptotic, post- measurement future (cf. 0]). By the 
equivalence principle, the analysis can be regarded as a 
model of gravity effects on entanglement. 

Our main results are: 

(i) In generic motion, particle creation in the moving 
cavity causes the entanglement to depend on time. 
This is in stark contrast to the previously-analysed 
relativistic situations (see [3, [f| for a small sample) 
where uniform acceleration is assumed to persist 
into the asymptotic future and the entanglement 
between inertial and accelerated observers is argued 
to be preserved in time. 

(ii) We give an analytic method for computing the en- 
tanglement for trajectories that consist of inertial 
and uniformly linearly accelerated segments in the 
small acceleration approximation. An advantage 



over the small amplitude approximations custom- 
ary with the dynamical Casimir effect Q is that 
the segment durations and the travelled distances 
are arbitrary. The entanglement remains constant 
within each segment, but we find that it does de- 
pend on the changes in the acceleration and on the 
time intervals between these changes. We present 
explicit results for three sample scenarios for a 
massless field in (1 + 1) dimensions, finding in par- 
ticular that in this approximation any degradation 
caused by the accelerated segments can be undone 
by fine-tuning the durations of the intermediate in- 
ertial segments. 

(iii) In (3+1) dimensions, the entanglement degradation 
has an observable magnitude for massless trans- 
verse quanta of optical wavelength at Earth gravity 
acceleration and for kaon mass quanta already at 
microgravity acceleration. The effect should hence 
be detectable in space-based experiments, where 
it would in particular test whether an equivalence 
principle between acceleration and a gravitational 
field holds also in the context of quantum informa- 
tion. 



II. CAVITY PROTOTYPE CONFIGURATION 

Let Alice and Rob ( "Relativistic Bob" Q ) be observers 
in (1 + l)-dimcnsional Minkowski spacctime, each car- 
rying a cavity that contains an uncharged scalar field 
of mass /i > with Dirichlet boundary conditions. Al- 
ice and Rob are initially inertial with vanishing relative 
velocity, and each cavity has length 6 > 0. The field 
modes in each cavity arc discrete, indexed by the quan- 
tum number n € {1,2,...} and having the frequencies 
uj n := y/AI 2 + TT 2 n 2 /S where M := /xS (we set c = h = 1). 

Let Alice and Rob initially prepare their two-cavity 
system in the maximally entangled Bell state \^f) = 
^ A \®)r + |l fe ) A \^k} R ), where the subscripts A and 
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Previously known as Fuentes-Guridi and Fuentcs-Schuller. 



R identify the cavity, |0) is the vacuum and |lfc) is the 
one-particle state with quantum number k. Experimen- 
tally, | might be prepared by allowing a single atom to 
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FIG. 1. Cavity trajectories in Minkowski space. Alice's cavity 
remains inertial. Rob's cavity is inertial in region I, acceler- 
ates in region II and is again inertial in region III. The figure 
shows the prototype case where Rob's acceleration is to the 
right and uniform throughout region II, and 77 is the duration 
of the acceleration in Rindler time ataxih(t / x) . 



emit an excitation of frequency u!k over a flight through 
the two cavities Q, and the assumption of a single k is 
experimentally justified if 5 is so small that cavity's fre- 
quency separation Lo n+ \ — lu„ is large compared with the 
frequency separations of the atom. 

We then allow Rob to undergo arbitrary acceleration 
for a finite interval of his proper time. After the acceler- 
ation Rob's cavity is again inertial and has proper length 
S in its new rest frame. Figure [1] shows the prototype 
case where Rob's acceleration is uniform while it lasts. 

Let U n , n = 1, 2, . . ., denote Rob's field modes that are 
of positive frequency iu n with respect to his proper time 
before the acceleration. Let U n , n = 1,2,..., similarly 
denote Rob's field modes that are of positive frequency 
oj n with respect to his proper time after the acceleration. 
The two sets of modes are related by the Bogoliubov 
transformation U m = ^2 n (a mn U n + (3 mn U*), where the 
star denotes complex conjugation and the coefficient ma- 
trices a and /? are determined by the motion of the cavity 
during the acceleration Q. Working in the Heisenberg 
picture, the state |\E') does not change in time, but for late 
time observations the early time states \0} R and |l/c)^j 
need to be expressed in terms of Rob's late time vacuum 
10)^ and the late time excitations on it, by formulas that 
involve the Bogoliubov coefficients |{|. In this sense, the 
acceleration creates particles in Rob's cavity. 

We regard the late time system as tripartite between 
Alice's cavity, the (late time) frequency u>k in Rob's cav- 
ity and the (late time) frequencies {uj n \ n ^ k} in Rob's 
cavity. As any excitations in the n 7^ k frequencies are 
entirely due to the acceleration, we regard these frequen- 
cies as the environment. We ask: Has the entanglement 
between Alice's cavity and the frequency Uk in Rob's cav- 
ity been degraded, from the maximal value it had before 
Rob's acceleration? 

We quantify the entanglement by the negativity, the 
widely- used entanglement monotone defined by Af(p) := 
— 2^ Ai<0 Ai, where the reduced density matrix p is ob- 



tained by tracing the full density matrix |\&) ($| over 
Rob's late time frequencies {w„ | n ^ k} and Ai are 
the eigenvalues of the partial transpose of p. Af has the 
advantage of being easy to compute in bipartite systems 
of arbitrary dimension [l(| • The closely-related logarith- 
mic negativity, E\a := ln(l +Af), is an upper bound on 
the distillable entanglement Ed and is operationally in- 
terpreted as the entanglement cost Eq under operations 
preserving the positivity of partial transpose [llj . In this 
respect, the entanglement quantification based on nega- 
tivity nicely interpolates between the two canonical (yet 
generically difficult to compute) extremal entanglement 
measures Ed and Eq [HJ ■ A sample of recent entangle- 
ment analyses utilising negativity can be found in [13J. 

The situation covering all the scenarios below is when 
a = diag(2i, z 2 , ■ ■ ■) + /ia (1) + 0(h 2 ) and /3 = + 
0(h 2 ), where h is a small parameter, the first-order coef- 
ficient matrices c^ 1 ' and are off-diagonal, and each 
z n has unit modulus. We find that the first correction to 
p occurs in order ft, 2 , the partial transpose is to this order 
an 8 x 8 matrix with exactly one negative eigenvalue, and 
the order h 2 formula for the negativity reads 



^ 2 E 



is (i) r 

\f J nk I 



(1) 



where the prime on the sum means that the term n = k 
is omitted. 



III. MASSLESS FIELD 

We now specialise to a massless field. Let I, II and 
III denote respectively the initial inertial region, the re- 
gion of acceleration and the final inertial region. As a 
first travel scenario, let the acceleration in region II be 
uniform, in the sense that Rob's cavity is dragged to the 
right by the boost Killing vector £ := xd t + td x as shown 
in Figure [TJ Let the proper acceleration at the centre of 
the cavity be h/5, where the dimensionless positive pa- 
rameter h must satisfy ft, < 2 in order that the proper 
acceleration at the left end of the cavity is finite. In re- 
gion II, the field modes that are positive frequency with 
respect to £ are then a discrete set V n with n S {1, 2, . . .} 
and their frequencies with respect to the proper time r at 
the centre of the cavity are O n = (vr/in)/[2<5atanh(/i/2)]. 

The Bogoliubov transformation from_f U n } to {V n } can 
be computed by standard techniques [8[ at the junction 
t = 0. The coefficient matrices, which we denote by a a 
and o/3, have small h expansions that begin 



1- ^Tr 2 n 2 h 2 + 0(h A ), 



(2a) 



(-1 + (-i) m ~ n ) 

— — '—^h + Oih 2 ) (m^n), 




7T 2 (m + n) 



h + 0(h 2 ). 



(2b) 
(2c) 



3 



The Bogoliubov transformation from region I to 
region III can now be written as the composi- 
tion of three individual transformations. The first 
is with the coefficient matrices ( Q a, a /3) from I 
to II. The second is with the coefficient matrices 
(diag(p,p 2 ,p 3 ,...),diag(p" 1 ,^ 2 ,p~ 3 ,...)), where p := 
cxp(if2if) and f is the duration of the acceleration in r: 
this undoes the phases that the modes V n develop over 
region II. The third is the inverse of the first, from II 
to III, with the coefficient matrices ( at,— Q/ 8 ). Col- 
lecting, we find from ([1} that the negativity M\ is given 
to order h 2 by 



the centre of the cavity between sending and recapturing 
a light ray that bounces off each wall once. 

As a second travel scenario, suppose that Rob blasts 
off as above, coasts inertially for proper time f ' and then 
performs a braking manoeuvre that is the reverse of the 
initial acceleration, bringing him to rest (at, say, Alpha 
Ccntauri). Composing the segments as above, we see that 
the negativity A/2 is periodic in f ' with period 25. Noting 
that for leftward acceleration holds with negative h, 
we find to order h 2 the formula 



A/- 2 = ±-/i 2 5> fer \p u 



-Am) 



l+2r 



M = \-h 2 Y j akr\ P 1+2r -l\ 2 

r=0 

= ±-2[Q(k,l)-Q(k :P )]h 2 , 



(3) 



where 



4/; 



Q( n , z ) ■= —r Re Li 6 (z) - — Li 6 (z 



1 



64 



_ 00 



+2r\ 



(l + 2r) 5 (l + 2r) 6 



(4) 



and the coefficients in the expansion Q(n,z) 



u-nr Re [z ) . Here Li@ is the polylogarithm [141 ] . 
the square brackets in the lower limit of the sum in Q 
denote the integer part, and a nr are all strictly positive. 
The sum term in ^ is 0(n~ 3 ) as n — > 00, and numerics 
shows that the sum term contribution to Q{n, z) is less 
than 1.1% already for n = 1 and less than 0.25% for 
n > 2. 

A/i © is periodic in f with period 27rf27 1 and attains 
its unique minimum at half-period. A plot is shown in 
Figure [2] The reason for the periodicity is that the full 
time evolution of the field in Rob's cavity during the ac- 
celerated segment is periodic with this period since the 
frequencies £l n are integer multiples of the fundamen- 
tal frequency fii. Ai is therefore periodic not just in 
the small h approximation of ([3]) but exactly for arbi- 
trary h. More generally, the same periodicity occurs for 
all cavity trajectories that contain a uniformly acceler- 
ated segment. We note that the period can be written as 
25 (h/2)~ 1 atanh(ft,/2): this is the proper time elapsed at 



= § - 2[2Q(fc, 1) - 2Q(k, P ) + Q(k,p') 
-2Q(k,pp') + Q(k,p 2 p , )]h 2 , 



(5) 



where p' :— exp(i7rr'/<5J. In addition to displaying the 
periodicities in f and f', ((5]) shows that the coefficient 
of h 2 vanishes iff p = 1 or pp' = 1 . This implies that to 
order h 2 , any entanglement degradation caused by the 
accelerated segments can be cancelled by fine-tuning the 
duration of the coasting segment. A plot is shown in 
Figure [3] 

As a third scenario, suppose Rob travels to Alpha Cen- 
tauri as above, rests there for proper time f" and then 
returns to Alice by reversing the outward manoeuvres. 
Again composing the segments, we see that the negativ- 
ity A3 is periodic in f " with period 25, and to order h 2 
we find 



A3 = I - h 2 ^ 1| V) 1 ^ - 1 



r=0 



X Wp'p") 



(6) 



where p" := exp(iwf"/6^, and the sum can be expressed 
as a sum of 14 Qs if desired [cf. the final expressions in 
((3|) and i[5])]. The periodicites in f , ?' and f" are manifest 
in ([6]). The coefficient of h 2 vanishes iff p = 1, pp' = 1 
or p 2 p'p" — 1, so that to order h 2 any entanglement 
degradation caused by the accelerated segments can be 
cancelled by fine-tuning the duration of either of the in- 
dependent inertial segments. Selected plots are shown in 
Figure |U 
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FIG. 2. The plot shows (| - Ai) h~ 2 with k = 1 as a function 
of it := Qif over the full period < u < 2tv. 



FIG. 3. The plot shows (| — A2) h 2 with k = 1 as a function 
of u := Clif and v := irf ' /8 over the full periods < u < 2ir 
and < v < 2n. Note the zeroes at u = mod 2-n and at 
u + v = mod 2-7T. 
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FIG. 4. The plots show (± - W 3 ) hT 2 with k = 1 as a function 
of u := Qif and v := nf'/S for < u < 2tt and < v < 2n 
when (from top to bottom) f" — 0, 25/3 and 4<5/3. 

Four comments arc in order. First, should one wish 
to consider noninertial initial or final states for Rob, our 
small h analysis is applicable whenever the assumptions 
leading to formula ([!]) still hold. For example, in a kick- 



and we have again verified that the Bogoliubov identities 
are satisfied pcrturbatively to order h 2 . The perturbative 
treatment is now valid provided h 1 and KM 2 < 100, 
allowing the possibility that M may be large. When 



start scenario that contains just regions I and II of Fig- 
ured! so that Rob's final state is uniformly accelerating, 
we find A/kick = | — Q(k, l)h 2 . 

Second, the validity of our perturbative treatment re- 
quires the negativity to remain close to its initial value \ , 
which in our scenarios happens when kh -C 1 . As the ex- 
pansions ([2]) arc not uniform in the indices, the treatment 
could potentially have missed even in this regime effects 
due to very high energy modes. However, we have veri- 
fied that when the h 2 terms are included in the expan- 
sions ([2]), these expansions satisfy the Bogoliubov iden- 
tities [8[ perturbatively to order h 2 and the products of 
the order h matrices in the identities are unconditionally 
convergent. This gives confidence in our order h 2 neg- 
ativity formulas, whose infinite sums come from similar 
products of order h matrices. 

Third, the matrices can be self-consistently trun- 
cated to the lowest 2x2 block provided the rows and 
columns are renormaliscd by suitable factors of the form 
1 + 0(h 2 ) to preserve the Bogoliubov identities to or- 
der h 2 . Taking Rob's initial excitation to be in the 
lower frequency, we find that all the above negativity re- 
sults hold with the replacement Q(l,z) — > aioRe(z) + 
iaii Re(z 3 ), and the error in this replacement is less 
than 0.7%. The high frequency effects on the entangle- 
ment are hence very strongly suppressed. 

Fourth, the analysis can be adapted to a fcrmionic field 
[lH and to scenarios where mode entanglement is gener- 
ated from an initially unentangled state fl6j . 



IV. MASSIVE FIELD 

For a massive field the frequencies are not uniformly 
spaced and the negativity is no longer periodic in the 
durations of the inertial and uniformly accelerated seg- 
ments. The massive counterparts of the expansions ([5]) 
can be found using uniform asymptotic expansions of 
modified Bessel functions [3, [lTj , with the result 



(7a) 
(7b) 
(7c) 

I 

k <C M, a qualitatively new feature is that the order 
h contribution in Q a is proportional to M 2 , resulting in 
an overall enhancement factor M 4 in the negativity. In 
the travel scenario with one accelerated segment, the neg- 



h + 0(h 2 ) (m^n), 
h + 0(h 2 ) (m#n), 



240 120 24QTr 2 n 2 96tt 4 ti 4 
2mn(~l + (-l) m " n ) [ir 2 {n 2 + 3m 2 ) + AM 2 } (M 2 + 7r 2 n 2 ) 1/4 

7T 4 (m 2 -n 2 f(M 2 +n 2 m 2 ) 1/4 
2mn(-l + (-l) m ~") [ir 2 {m 2 + in 2 ) + AM 2 ] (M 2 + 7r 2 m 2 ) 1/4 
7r 4 (m 2 -n 2 ) 3 (M 2 + 7r 2 n 2 ) 1/4 
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FIG. 5. The plots show (± - Mi) h~ 2 M- 4 © for M = 10 3 as 
a function of nf /(AM 8), in the upper figure with k = 1, 2, 3, 4 
(solid, dashed, dash-dotted, wide-dashed) and in the lower 
figure with k = 30. 



ativity takes in this limit the form 



, ,,„, 256k 2 



i 



(k 2 -n 2 y 
(VM 2 + 7r 2 fc 2 - yjM 2 + n 2 n 2 ) (f/S) } 



(8) 



where the double prime means that the sum is over pos- 
itive n with n = k + 1 mod 2. A/i <j8j> is approximately 
periodic in f with period iMS/ft, but it contains also sig- 
nificant higher frequency components. Plots are shown 
in Figure [5] 



cavity of length 8 = 10 m and acceleration 10ms -2 , an 
effect of observable magnitude can be achieved by trap- 
ping quanta of optical wavelengths provided the momen- 
tum is highly transverse so that k -C M « 10 8 . Were 
it possible to trap and stabilise massive quanta of kaon 
mass fi = 10~ 27 kg in a cavity of length 8 = 10 cm, the 
effect would become observable already at the extreme 
microgravity acceleration of 10 _10 ms -2 . 

These estimates suggest that experimental verification 
of the effect is feasible, but they also suggest that ac- 
celerations of gravitational origin should be properly ac- 
counted for, both in the theoretical analysis and in the 
experimental setup. On the theoretical side, we antici- 
pate that the core properties of our analysis extend to 
nonlinear acceleration and to weakly curved spacctime. 
To achieve experimental control, especially on accelera- 
tions on gravitational origin, the experiment may need 
to be performed by entangling cavities in spaceships. An 
experiment within a single spaceship could test the rela- 
tive acceleration effect analysed in this paper. An experi- 
ment using cavities in two and possibly widely separated 
spaceships would test whether an equivalence principle 
between acceleration and a gravitational field holds also 
in the context of quantum information. An experimental 
confirmation that gravity degrades entanglement would 
indeed provide a novel addition to the currently scarce 
experimental evidence on quantum phenomena involving 
gravity [l8j . 
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tive (1 + l)-dimensional mass. For a massless field in a 
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